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2 ' Abstract 

'NT 

. We compute the ground state energy of a massive scalar field in the 

[ background of a cylindrical shell whose potential is given by a delta func- 

tion. The zero point energy is expressed in terms of the Jost function of 
the related scattering problem, the renormalization is performed with the 
help of the heat kernel expansion. The energy is found to be negative for 
^ ' attractive and for repulsive backgrounds as well. 

> 

X 1 Introduction 

The response of vacuum to classical external fields or constraints has been studied 
since more than fifty years. The manifestation of this response are macroscopic 
forces exerted on boundaries or surfaces where the potential is concentrated. The 
most famous example is the Casimir effect IQ , which has been recently measured 
in laboratory with a precision of nearly 1% Notwithstanding the large size of 
theoretical work and the experimental efforts many question remain still open. 
We have in fact by this far no exact knowledge of the law governing the sign of 
the Casimir energy, how exactly the features of the boundary and its geometry 
influence the response of vacuum. An investigation on the cylindrical geometry 
has begun in after the suggestion that a cylinder, as a kind of intermediate 
shape between the parallel plates and the sphere, could possess a zero Casimir 
stress 0. However paper |Q showed that a perfectly conducting cylinder in the 
electromagnetic vacuum has a negative Casimir energy. A study of the dielectric 
cylinder has been recently carried out in |^ and in with the interesting result 
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of a vanishing vacuum energy in the dilute case. Some hght could be further 
shed into the problem with the study of other quantum fields and other types of 
background potentials. In this paper we investigate a cylindrical shell with radius 
R, having a delta function 6{r — i?) as a potential profile. A delta function is 
an idealization of a strongly concentrated potential. Unlike the case of Dirichlet 
boundary conditions the field is continuous on the surface of the shell, therefore 
the delta function represents a model for a boundary which is not completely 
hard, but which becomes transparent at high frequencies. This model can be 
considered as the "scalar" version of the dielectric background. A non singular 
potential would be, of course, more realistic under the physical point of view, 
but the calculations would be considerably more complicated. A calculation with 
the delta function has been already carried out for a spherical shell [0 and for 
a magnetic flux tube in the vacuum of spinor and scalar fields [Q. We hope 
that our work can contribute in the understanding of vacuum as a fundamental 
aspect of quantum field theory. The cylindrical shell could, besides, find physical 
applications in the calculation of the quantum corrections to vortices in QCD or in 
the electroweak theory. Another interesting perspective is in the recent discovery 
of the so called nanotubes §], fl^, which are large carbon molecules generated 
in laboratory offering the intriguing possibility of measuring quantum effects on 
small cylindrical objects. It is also of interest the comparison of the problem of 
the cylindrical shell with the analogous problem of the magnetic flux tube ||^. 
The organization of this paper is as follow: in the next section we collect the 
basic formulas necessary for the calculation and we present the renormalization 
procedure, in section 3 and 4 the major part of the analytical calculation is 
performed. In section 5 the numerical results are given and in section 6 these 
results are briefly discussed. 

2 Vacuum energy in terms of the Jost function 
and renormalization 

Let us work with a real massive scalar field t) with mass m and let us quantize 
it in the background of a cylindrical potential. The field equation in cylindrical 
coordinates r, (f, z, after separation of the variables, reads 

(^pl-m^-pl-^--Vir) + ^dr + df^MPo,Pz,r) =0, (1) 

where is the momentum four vector being pz its component along the longi- 
tudinal axis of the cylinder and I is the angular momentum quantum number. 
V{r) is the background potential given by 

V{r) = ^6{r - R) , (2) 
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it represents an infinitely thin cylindrical shell whose profile is a delta function. 
The shell has a circular section of radius R and it extends from —oo to +00 in the 
z direction, a is a dimensionless parameter giving the strength of the potential. 

Considering eq.(l) as a scattering problem, we choose the "regular solution" 
which is given by 

0(r) = JKA;r)e(i?-r) + ^(/KA;)i7P(A;r) + /;(A;)i7/^)(A;r))0(r-i?) , (3) 

where k = a/Pq ~ ~ pIj Jm{kr) is a Bessel function of the first kind, H[^\kr) 
and Hf\kr) are the Hankel functions of the first and of the second kind and 
the coefficients fi{k) and fl{k) are a Jost function and its complex conjugate 
respectively. Q{R — r) is a theta function. The field is therefore free in the 
regions < r < i? and R < r < 00. At r = R the field is continuous and we have 
the following matching conditions 

0'(i? + O)-0'(i?-O)= 

(PiR + 0) = - 0) , ^ ^ 

where the prime indicates derivative with respect to r. 

The quantum field is in its vacuum state, its energy is given by half of the sum 
over all possible eigenvalues u!(^n) of the Hamilton operator related with the wave 
equation (1). We define a regularized vacuum energy 

(n) 

where s is the regularization parameter which we will put to zero after the renor- 
malization, /x is a mass parameter introduced for dimensional reasons and (n) 
includes all possible quantum numbers. We will calculate the energy density per 
unit length of the cylinder given by 

°° (n) 

where e(„) are the eigenvalues of the operator contained in (1) without Pz and 
with k = \/p\— va^ ■ We carry out the integration over and we arrive at 

(n) 

Following a well known procedure]] we transform the sum in (7) into a contour 
integral and, dropping the Minkowski space contribution, we find 

e = -jCs / dk{k^ + mY~'d,lnUik) , (8) 



^ The procedure is explained with details in [ pl| . 
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where = (1 + s(— 1 + 21n(2/i)))/(27r) and fi{ik) is the Jost function defined 
in (3) on the imaginary axis. It contains all the information about the back- 
ground potential under examination. We will find fi{ik) explicitly in the follow- 
ing section. The energy defined in (8) is renormalized by direct subtraction of its 
divergent part 

^ren ^ ^div (9) 

with the normalization condition demanding that the vacuum fiuctuations vanish 
for a field of infinite mass 

lim Sren = . (10) 



m— >oo 



This normalization condition fixes the arbitrariness of the mass parameter /i (for 
more details see |12]. The condition does not apply to massless fields). The 
isolation of £div is done with the use the heat-kernel expansion of the ground 
state energy 

S = E ^^-^^P^m^-^^^-^^A, , , = 0, 1 1, ... , (11) 
i 

where the Aj are the heat kernel coefficients related to the background. Then we 
define 

1 \ Tn? 

Tfl 







+ In 


4^2 


647r2 




IV? 






+ In 


4/i2 


' 327r2 




711? 


- 1 




+ In 


4^2 


327r2 




IT? 



- 2J ^2 , (12) 

which includes all the pole terms and all the terms proportional to non-negative 
powers of the mass. 

In order to perform the analytical continuation s — > we split the renormal- 
ized vacuum energy (9) into two parts 



with 



and 



^/ = -i^^Ey dk[k'-m']'-^ — [\nm)- In frm] (14) 



/=— oo ' 



1 °° d 
£as = —.CsY^ / dk[k''-m^Y-'—\nfr{ik)-£d^.. 
4 ,.Jm ok 



(15) 



Z=— oo 
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Here fj^^{ik) is a portion of the uniform asymptotic expansion of the Jost function 
which must include as many orders in / as it is necessary to have 



(16) 



in the hmit / — oo, — oo with fixed l/k, which is sufficient to let the sum 
and the integral in (14) converge at s = 0. The splitting proposed in (13) leaves 
the quantity Sren unchanged, while it permits the substitution s = in the finite 
part Sf. 

3 The Jost function and its asymptotics 

We insert solution (3) into (4) and we find 



MkR) = i mHl'\kR) + ft{k)H\"{kR) 



(1), 



fl{k)Hf\kR) + fl{k)H\'\kR) ) \r=R = fJlikR) + idrJlikr)) \r=R . 

(17) 

The system can be solved for fi{k), with the help of the Wronskian determinant 



r(l)/ 



of the Hankel functions [|13[ , the result is 

fi{k) = -\i'KkR(ji+r{kR)H^^\kR)-Ji{kR)H\'^^{kR) 



r(i) 



a 



-^J,{kR)Hi'\kR) 
The corresponding Jost function on the imaginary axis can be written in terms 



of the modified Bessel I and K functions, again with the help of |13 



Mzk) = l + aIi{kR)Ki{kR). 



(19) 



From the Jost function (19) one arrives at the uniform asymptotic expansions 
fi''^{ik) for positive I and f^^~{ik) for negative I and at the asymptotic expansion 
fQ^{ik) for / = 0, with the help of the asymptotics of the Bessel I and K functions 

Since the asymptotic Jost 



for large indices and large arguments available on 113 



function consists of three different contributions, the sum over / appearing in (14) 
and (15) must also be distinguished in three contributions: a sum over negative 
/, a sum over positive / and a contribution coming from / = 0. The first two 
contributions can be summed up analytically in the following way 



...ln/r(^A:)... + X;-ln/rW- = J] -(1^ /r"- W + ln / 



l=—oo 



1=1 



{ik))... , 

(20) 
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where the dots represent, for simphcity, the rest of the functions in (14) and (15). 
Then, eq.(14) and (15) are rewritten in the form 



--Csj dk[e-m^]'-^— In fS'{zk)-S, 



(21) 



and 



-1 CO 

Sf = dk[e-m']—{2lnm)- lnfrH^k)) 

-—j dk[e-m']— (In foitk)- In fS'itk)) , (22) 

where In f°-^^{ik) = In (iA;) + In f^^i~{ik) and we have used the property 
fi{ik) = f^i{ik) of eq.(19). 

Taking the logarithm of the uniform asymptotics of the modified Bessel func- 
tions and re-expanding in negative powers of the variable / (see for details on 
this procedure, see also appendix A) we find 

^ 2kR ~ Wb^ ' ^^^^ 
Inr^itk) = Y.Y.^n~, (24) 

n=l j 

where t = (1 + {kR)/lY)^ and the non-vanishing coefficients are 

Xi,i = a , X2,2 = -aV4 , 

^3,3 = a/8 + ay 12 , ^3,5 = -3a/4 , (25) 
X3J = 5a/8 . 

In this definition we have included 3 orders in I which are sufficient to satisfy 
condition (16). Substituting (23) and (24) in (21) and (22) we find 

1 r°° f) / _ ^ 



1=1 \n=l j 



P 



and 



-ic. /; ,kie - - ^) - (26) 

^/ = -^E/ ^M^^--^]^ 21n/K.A:)-$:$:X.,- 

Z=l \ n=\ j / 



We call the first addend in (27) Efi and the second £^/o, that is Ef — £fi + £fo. 



4 Asymptotic part of the energy 

We go forward with an analytical simplification of (26). We call the second 
addend in (26) Saso, it can be immediately calculated, giving 



£asO — 



am 



a 



647ri?2 V s 



+ In 



(28) 



A simphfication of the first addend in (26) can be achieved with the Abel-Plana 
formula (appendix B) which transforms the sum over I into an integral. Then, 
the first addend in (26) turns out to be the sum of three contributions 



S, 



asl 



am 

Ten 

a^m 
^32R 



l+ln(V 

s \ m'^ 



£as2 — 



am 



+ 



a 



SttR 647ri?2 Vs 



+ In 



m^ 



+ 



a 



- + In 
s 



a 



a 



647rmi?3 967rmR^ 



(29) 
(30) 



S. 



as 3 



a 

hAmR) — 



a' 



32i?2 

HmR) + ^ himR) , 



'^^^^)+(l6^ + 2^) 



(31) 



where the functions hi{x) are convergent integrals over / given explicitly in ap- 
pendix B. In appendix B the reader can also find more details about the derivation 
of the three contributions (29), (30) and (31). £asi and Eas2 contain all the pole 
terms (all the divergences of the vacuum energy) plus the terms proportional to 
non-negative powers of the mass (which do not satisfy the normalization condi- 
tion). All these terms are subtracted and are used to calculate the heat-kernel 
coefficients by means of definition (12) for Sdiv Below we give the heat kernel 
coefficients which we calculated up to the coefficient (adding four more orders 
in In f^^^{ik)), in the hope that they will be of use for future investigations on 
the same background 



As 
A, 



= 



Ai/2 = 



— 27rQ; 

210R4 
(64Q3+52aS+39a^)7r 
1638aR^ 



i3/2 
h/2 



2R 

_(3af+4aV!^ 

(81a^+120a"^+128a'^)7r3/^ 
24576i?s 



(32) 
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One can note how all the integer heat-kernel coefficients depend on odd powers of 
the coupling constant, while the half- integer coefficients depend on even powers 
of a. The same feature is present in the heat kernel coefficients of a 5-potential 
spherical shell 0. 

We perform the subtraction of the divergent portion and we obtain the final 
result 

S-riH = ^ hi(mR) — hoimR) 



8nR^ /^4(mi?) + ^ h,{mR) - - . (33) 

5 Finite part of the energy and numerical re- 
sults 

The quantity (27) cannot be analytically simplified. We integrate Sfi and Sfo by 
parts and make the substitution k ^ k/R to get an exphcit dependence on R 



1=1 \ \n=l j y 

Sfo = ^l/kk ^ln/o(zfc)k.^./^-- + —J 



(34) 



For small values of R, Sji behaves like R ^, while Sjo has a logarithmic behaviour 

ln(mi?) 

32^ ■ ^^'^ 

We note that in the limit a ^ 0, the logarithm of the Jost function fi{ik) can be 
expanded in powers of a, giving (from eq.(19)) 

\n{l + aIi{kR)Ki{kR)) aIi{kR)Ki{kR) + ^a^I^{kR)Kf{kR) + 0{a^), 

then, the summation over / of the leading term of this expansion could be analyt- 
ically performed following a method recently proposed in P] . This would give us 
a ffist order approximation of Sfi. However we will restrict us here to a numerical 
calculation of the sum in (34), which can be performed for small as well as for 
large values of a. 

To find the asymptotic behaviour of Sas we rewrite eq.(33) in the following 
form ^ 

Sas = 2^^2 + a^W2{mR) + a^wzimR)] , (36) 
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where the functions wi{x), W2{x) and w^i^x) are given by 



111 1 

wi{x) = { -hi{x) + -hs{x) - -h4{x) + —h^ix) - 



W2{x) = y- — h2{x] 



-3(.) = (^ft3(.) - ^) . (37) 

As mentioned above, the functions hi{x) are quickly converging integrals. The 
behaviour of Sas is governed by the functions ^1,2,3(2^) and by the value of a. For 
— > we find 

wi{x) ~ 0.0000868 + 0{x) , 

W2{x) ~ —\nx + 0.115 + 0{x) , 
16 

W3{x) ~ ^Inx + 0.00479 + 0{x) . (38) 

Thus, for small values of the radius, Sas is proportional to Ini?. Summing the 
contribution coming from (35) and (38), we find that the leading term of the 
renormalized energy, for R 0, is 

\n{mR) 

^ A8nR- ■ ^^^^ 



This result is in agreement with the prediction (section II of paper jT^ ) 



lim £ren ~ HmR) . (40) 

In the limit R ^ 00 the behaviour of the renormalized energy is determined by 
the first non-vanishing heat kernel coefficient after A2. From (32) and eq.(ll) we 
arrive at 

/ 1 \ 

lim Sren ~ - " - , , J + ■ (41) 



R^oo 2048mi?3 1536mi?3 

We have numerically evaluated the quantities Sas, Sfu ^fo and Sren as func- 
tions of R, fixing the values of the mass to 1. It turned out to be necessary 
to sum 20 terms in the variable / and to integrate up to 1000 in the variable k 
to obtain "stable" numerical values for the energy. Below we give the plots of 
the various contributions to the vacuum energy and the complete renormalized 
energy for different values of the potential strength. For a < we found the 
renormalized vacuum energy to posses a small imaginary part becoming larger 
when a approaches —oo. Particle creation accounts for this contribution. It 
starts when the attractive potential of the shell becomes over-critical, that is 
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when e < —m < 0, where e is the energy of the bound state; in this case the 
effective action of the system acquires an imaginary part0, however a detailed 
discussion of this aspect of the theory is beyond the purpose of the present pa- 
per. It should only be said that in the plots traced for negative values of a (Fig. 
4 and 5) the energy is to be intended as real part of. 



6 Conclusions 

We calculated the vacuum energy of a scalar field in the background of a cylindri- 
cal semi-transparent shell. The formulas for the energy density per unit length of 
the shell are given by equations (33) and (34). The heat kernel coefficients (32) 
for a cylindrical potential containing a delta function, are also a relevant part of 
the results. A discussion of the sign of the vacuum energy is possible by means of 
the data given in the numerical section of this paper. The energy is found to be 
negative in the background of repulsive potentials {a > 0) for every finite value 
of the radius (Fig. 3). This conclusion sets a closeness between the model studied 
here and that of a conducting cylinder [Q and of a dielectric cylinder |^ in the 
electromagnetic vacuum. The latter models possess also a negative Casimir en- 
ergy. There is also a resemblance with the 5-potential spherical shell investigated 
in paper [0, where a negative energy was observed for large repulsive potentials. 

In the background of attractive potentials {a < 0) the energy is negative 
on almost all the i?-axis. As we see in Fig. 5, the energy becomes positive for 
very small values of the radius. This happens when the logarithmic term in (39) 
compensates for the term proportional to — in (41) which is dominant for 
small values of a in the region R > 1. The thin region of positive asymptotically 
increasing energy appearing in Fig. 5 is reasonably out of physical applicability. 
For both repulsive and attractive potentials the renormalized vacuum energy of 
the semi-transparent cylinder goes logarithmically to ±oo in the limit i? — 0. 
This logarithmic behaviour was also found in the semi-transparent spherical shell 
0, however it was not found in the ^-potential fiux tube P|, which has many 
feature in common with the background investigated here. This can be explained 
by means of the heat-kernel coefficient A2, which is a non zero coefficient here 
as well as in 0. It is also interesting to note how in expression (32) for the 
coefficient A2 the contributions proportional to a and have cancelled and 
only a term proportional to is present. The cancellation of the lower powers 



of the potential strength was also observed in 0] and in paper |T^, where a 
theta function profile in a dielectric spherical shell is examined. It confirms the 
observation that singular profiles show weaker divergences than smoother, less 
singular profiles. 



^For the theoretical foundations of this phenomenon see jl^l- In the context of vacuum 
energy the appearance of an imaginary part in the rcnormahzcd energy was observed for instance 
in M. 
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8 Appendix A: Asymptotic expansions of the 
modified Bessel functions 

The expansion in negative powers of the parameter / (eq. (24)) is obtained from 
the following expansions 

1 ^ 

h+am) ~ ^=cxp{ V r"5,(^,«,^)} , (42) 
Ki+aikR) ~ y|^exp{ J] r"5;,(n,a,i)} (43) 

* n=-l 

where a takes the values 1 and and the functions Si{n, a, t) and Sxin, a, t) up 
to the third order are given by 



r< / 1 1 +t 

Si{Q,a,t) = -mt — -In 



2 2 \l-t 

Si(l,a,t) = -^{-3 + 12a'^ + 12at + 5f) , 

Si{2, a, t) = —[8aH + 12a^(-l + 2t^) + a{-26t + 30t^) + 3(1 - 6t^ + 5t^)] , 
48 

,5/(3, a, t) = — (((25 - 104a^ + 16a^)t^)/3 + 16a(-7 + 4a2)i^ 
128 

-(531/5 - 224a2 + 16a'^)f^ - (32a(-33 + 8a^)t^)/3 

-(-221 + 200a2)t^ - 240at^ - (1105t^)/9) ; (44) 



^ / N 1 , a , /I + t 
5^(0,0,^) = -lni + -ln' 



2 2 Vl-V ' 

5K(l,a,t) = -^{-3 + 12a^ + 12at + 5t^) , 

+2 

,5^(2, a, t) = :7^[-8a^^ + 12a^(-l + 2t^) + a{-26t + 30t^) + 3(1 - + 5t^)] , 
48 
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SK{3,a,t) = - — ((25 - 104a^ + 16a^)i^)/3 + 16a(-7 + 4a2)i^ 
128 

-(-531/5 + 224a^ - 16a^)t^ - (32a(-33 + Sa^)f)/3 

-(221 - 2Q0a^)t^ - 240at^ + (1105t^) /9 , (45) 

with i = (1 + (A;i?/0^)-i 



9 Appendix B: Calculation of the sum and of 
the integrals 

The transformation of the summation over I in (26) into an integral is done with 
the Abel-Plana formula 



In our case the function F{1) is 



/■oo / 3 

F{1)= dk{k^ + m^Y-'dAY,Y.^- 

-'^ \n=l i 



(47) 



To the first addend of cq.(46) corresponds the contribution Sasi to the second 
addend the contribution Sas2 ant to the third the contribution Sass- In ^asi the 
integrations over k and I are performed by means of the following formula 

giving directly result (29). For the integration over k in Sas2 and Sas3 the following 
formula has been used 



J m 



' * ie - ^Y-3jL . r(2-.)r(.^|-i)P^ _ ^^^^ 



which, applied to £as2-i leads directly to result (30). In Eas^ the integrand is more 
complicated and the integrals over / cannot be analytically calculated. Result (31) 
represents the maximum possible simplification, which we reached after several 
partial integrations. The functions hi{x) are 

dl 



12 



1 V' 



X 



1 - e2^' / 
1 1 



h{x) = I dl{ ) VP^x^ 



oo 



I yv' 



K{x) = I di{{^—^ \ -) 



h{x) = I dl{{{ ^—^ - - I VW^^ . (50) 
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Figure 3: Repulsive potential. The complete renormalized vacuum energy Sren{R) 
multiplied by ■ o;"^, for a = 0.3, a = 2.1 and a = 10. 
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Figure 5: Attractive potential. The complete renormalized vacuum energy 
£ren{R) multiplied by R"^ ■ a"'^, for a — —0.3. The peak beneath the lowest 
point of the curve is due to the presence of a small imaginary part in the renor- 
malized energy causing the function to be no more analytical. The energy starts 
to be complex at i? ~ 0.04. 
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